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A dimensionally reduced expression for the QCD fermion determinant at finite temperature and 
chemical potential is derived which sheds light on the determinant's dependence on these quan- 
tities. This is done via a partial zeta-regularisation, formally applying a general formula for the 
zeta-determinant of a differential operator in one variable with operator-valued coefficients. The 
resulting expression generalises the known one for the free fermion determinant, obtained via Mat- 
■ subara frequency summation, to the case of general background gauge field; moreover there is no 

' undetermined overall factor. Rigorous versions of the result are obtained in a continuous time — 

, lattice space setting. The determinant expression reduces to a remarkably simple form in the low 

' temperature limit. A program for how to use this to obtain insight into the QCD phase transition 

at zero temperature and nonzero density is outlined. 

<t^ . 11.15.Tk, 12. 38. Aw, 11. 15. Ha 

, QCD at finite temperature and density exhibits interesting phase structure; in particular, at sufficiently high 
temperature and/or density there is a transition from the usual confined, chiral symmetry-broken hadronic phase to 
, a deconfined, chirally symmetric phase where the quarks and gluons are "liberated" and form a quark-gluon plasma 
K*" (QGP). The Universe was very likely a QGP for a brief moment after the Big Bang, and such a phase may also exist 
at present in the cores of very dense stars (neutron/quark stars). Furthermore, it is a major aim of current heavy 
' ion collision experiments at RHIC and CERN to create a QGP for brief instants through the high-energy collisions of 
I , heavy nuclei. QCD at finite temperature and density is therefore an exciting field of major current research interest 
' [1]. The equilibrium properties are governed by the QCD grand-canonical partition function (referred to hereafter 
simply as the 'partition function'). A crucial quantity contained in this is the fermion determinant, which encodes 
' the dynamical fermion effects. The central problem in QCD thermodynamics is to determine the dependence of the 
I partition function on the temperature T and (quark) chemical potential fj, , and to this end it is highly desirable to 
' get as much information as possible on how the fermion determinant depends on these quantities. In this paper we 
Oh, derive a "dimensionally reduced" expression for the fermion determinant which sheds light on these dependencies. 
^ ■ The importance of the fermion determinant in this context has been demonstrated in Lattice QCD studies at 
'T"! I low temperature and nonzero chemical potential. It was found that when the fermion determinant is discarded (i.e 
K>" ■ in the "quenched approximation") the onset of the chiral phase transition in the quark number density and chiral 
condensate occurs at fio = ^tt/S (where is the lightest meson mass, i.e. the pion mass) rather than at the 
5^ ■ expected critical value /Zc = mj^/i (where is the lightest baryon mass, i.e. the nucleon mass) — see [2] and 
\ the references therein. The necessity of including the fermion determinant has also been demonstrated in another 
approach based on Random Matrix Theory [3]. The phase transition at zero temperature and nonzero density is of 
central interest since, e.g., it is expected to occur in the formation of sufficiently dense stars. Unfortunately, this 
region of the phase diagram is not currently accessible by the currently developed (non-quenched) lattice methods 
or perturbative methods [1]. However, there is an indication that it may be possible to obtain insight into this 
region by nonperturbative analytic techniques. Recently, T. Cohen has derived an interesting expression for the 
ratio detD A{y) / detD a{^) where detDA{^i) is the finite temperature fermion determinant in a background gauge field 
A and at chemical potential fj, [4]. In the low temperature limit it reduces to a remarkably simple form and the 
fermion determinant is seen to undergo a transition from /^-independence to /t-dependence at precisely the critical 
value Ho mentioned above [4]. The zero-temperature QCD partition function itself should undergo a transition at 
the above-mentioned critical value fic- Its /^-independence in the region < fJ- < fJ-c must therefore involve a 
subtle cancellation in the functional integral of e"'^"^^^'^'^) Ylfj^idetDAifJ- ,rnj) over gauge fields (where Sym{A-) is 
the Yang-Mills action, Nf the number of quark flavours and rrij their masses). This is intimately connected with 
the nature of the strong interactions, and in particular, confinement, in usual QCD (i.e. at T~fi = 0), since in the 
absence of these interactions we would simply have m7r/2 = m^v/S =lightest quark mass, i.e. = /tc- Therefore, 
understanding how this cancellation works is not only key to understanding the hadronic-to-QGP phase transition in 
QCD at zero-temperature and nonzero density but should also give insight into confinement in usual QCD. 
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Understanding this cancellation is therefore an interesting and important problem (dubbed the 'baryon Silver Blaze' 
problem in [4]). To make progress on this, a necessary first step is to supplement the determinant ratio expression of 
[4] with an expression for detDAin) itself, since we need to know its full gauge field-dependence in order to study the 
aforementioned cancellation in the functional integral over gauge fields. (A knowledge of detDAOj.)/d.etDA{0) does 
not suffice for this since we cannot a priori rule out the possiblity that there are important gauge field-dependent 
factors which drop out in the determinant ratio.) This provides a concrete motivation for the considerations in the 
present paper. The "dimensionally reduced" expression for detDAifi) that we obtain, (9) below, also reduces to a 
remarkably simple form in the low T limit ((17)-(18) below) and we identify a gauge field-dependent factor which 
drops out in the determinant ratio. Our results are derived both at the formal continuum level and rigorously in 
a continuous time — lattice space setting. (In the latter setting the evaluation requires certain choices to be made, 
resulting in different overall factors in the determinant expression, and it is a non-trivial question whether these are 
physically equivalent. This issue is intimately connected with a "universality anomaly" in Lattice QCD uncovered 
recently in [5].) The natural next step in this program, which is currently imder investigation, is to consider the strong 
coupling limit of the functional integral expression for the QCD partition function, in the low temperature limit with 
the fermion determinant expression obtained here, and compare it with the expression obtained in the strong coupling 
lattice Hamiltonian framework at finite chemical potential. (Note that this comparison makes sense in the continuous 
time — lattice space setting.) The latter has been studied, e.g., in [6] and the phase transition at the appropriate 
yttc has been explicitly demonstrated in that framework. Using the equivalence between the functional integral and 
Hamiltonian frameworks it should be possible to see how the aforementioned cancellation in the functional integral 
over gauge fields for Ho < jJL < comes about in the strong coupling limit. The final, and most challenging, step in 
the program will then be to extend the understanding of this cancellation to the case of general couplings. 

The "dimensional reduction" method introduced here is also potentially useful in other contexts and we briefly 
discuss several of these at the end of the paper. 

The QCD partition function with Nf quark flavours with masses mj and quark chemical potential can be expressed 
as a functional integral, 

N f 

J j=i J 

where (3 = 1/T. The fermion fields ipj , 'ipj (resp. gauge fields A^) are required to satisfy anti-periodic (resp. periodic) 
boundary conditions in the Euclidean "time" variable r 6 [0,/^]. 3- space is taken to be a box (whose volume V 
is taken to infinity in the final step of calculating physical quantities). The spacial b.c.'s will not play a role in 
our considerations and we leave them unspecified. The integrals over the fermion fields in (1) give the fermion 
determinants, i.e. for mj=m the integral gives detDA{l-t) where 

Da{ij.) = J^{d^ + A„) + m - i^'y4 (2) 

In the free field {A = 0) case, where detDo{iJ,) is the partition function of a free Dirac fermion gas, the dependence of 
the determinant on (3 and /i is well-known: An application of the standard Matsubara frequency summation method 
gives [7] 

detDoin) = Co n i^"" (1 + e-''(^-''))(l + e->'^'^+^^)f (3) 
p 

with ±E = ±-\/p^ + m^ being the (2-fold degenerate) energy eigenvalues of Hq = ■y4{jkdk + m) (summation 
over k = 1,2,3 is implied). From this, expressions for physical quantities such as the energy density {E)/V = 
— y-^ logdetDo |/3^=const and particle number density {Q)/V = yj§ji logdetDo are obtained; in particular one finds 
for massless fermion the well-known results {E)/V ^ ji^ and {Q)/V ~ /i^ in the large 13 limit. The Matsubara 
summation produces an undetermined overall factor Co in (3); this is inconsequential though since it does not involve 
(3 , jjb ov the energies -izE. 

We now consider the problem of generalising (3) to the case of arbitrary background gauge field. The Matsubara 
summation method can be used to obtain an expression for the ratio detD Aifj) / detD a{Q) [4]; however, it is of limited 
use for detDAifJ.) itself since, e.g., it produces an undetermined overall factor and we cannot exclude a priori the 
possibility that this factor may depend on the gauge field. Therefore we take a different approach. Regarding the 
spinor fields V'(x, r) as functions ^'(t) taking values in the vectorspace W = {■^(x)} of spinor fields living only in the 
spacial volume (i.e. the usual Hilbert space of quantum mechanical wave functions), we re-express Da{h) as 
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where Ha{t) :W is the Unear map defined by 

Ha{t)^{x) = [-yii'jkidk + Akix, r)) + m) + A^ix, t)]V(x) (5) 

In this way DA^fJ.) can be viewed as a differential operator in one variable r G [0, /3] acting on M^-valued anti-periodic 
functions. An expression for its determinant can then be obtained by application of a zeta-regularised determinant 
formula in [8] (see also [9] for related and overlapping results). Recall that the zeta-determinant [10] of an operator 
D is defined by det^D = "HjAi" = e-^^^^") where {A^} are the eigenvalues of D and Cd{s) is the zeta function 
defined by Cd{s) = X^- jf = Tr{D~^). Under certain ellipticity conditions Cois) is a well-defined smooth function 
of the complex parameter s wIicti Re{s) is sufficiently large, and can be analytically continued to a mcromorphic 
function in the whole complex plane which is regular at s = 0, so that Coi^) ^^'^ hence det^D are well-defined. The 
zeta-determinant formula of [8], in its most straightforward form, is for elliptic differential operators in one variable 
t G [0, 1] acting on periodic functions taking values in some vectorspace. To apply it in the present case, we note that 
the spectrum of Da{ij) is the same as that of the operator Da(p) = 74 ( + HAiPt) — /i — ^) acting on W^-valued 
periodic functions '^{t) , t G [0, 1] : It is easy to check that \E'(r) is an eigenfunction for Da{ij) with eigenvalue A if 
and only if 4i{t) = e"^* is an eigenfunction for DAin) with the same eigenvalue. Hence detDA{fi) = detDA{n), 

and the result of [8] can be applied to evaluate the latter determinant. Writing Da{ij) = 7^ + LQ{t) where 
Lii{t) = ^74 , Lo{t) = 74 {HA{l3t) — /U — ^) , and formally setting N = dimW, Theorem 1 of [8] gives 

detDAilA = deWAilA = (-1)"^ Se{Li,Lo) R{Lu io) detil - U{Lx, Lo)) (6) 

where the ingredients are as follows. Consider the equation DA{n)'^{t) = without boundary conditions on 
^(f). The solutions are determined from their initial values via an evolution operator: ^{t) = U{t)'i{0). Then 
U{L\,Lq) := U{1) = Te~'Jo Lo{t)dt ^y=^.ordering). The remaining ingredients are given by R{L\,Lq) := 

^ir,Tr(L,{t)-'Lo(t))dt Se{Li,Lo) := (deire)e'./"o ^'■(^^W^^W^^W)'^* with Teit) as defined in the following. 

Note that for A G C the definition of A~* depends on a choice of cut in the complex plane; the same is true for Cd{s), 
and hence the zeta-determinant det^D also depends on such a choice. In the formula (6) this choice is represented 
by 6. Specifically, choose ^ G R so that the line {re*^ | r G R} does not contain any eigenvalue of Li{t). Then the 
argument a of an eigenvalue of Li{i) satisfies cither 9<a<9 + TTOT6-'r'rr<a<9 + 27r. Let n^(t) and n^(t) 
denote the projections onto the subspaces of W spanned by the eigenvectors of Li{t) for which the arguments of the 
corresponding eigenvalues are in the intervals ]9,9 + 'k[ and ]9 + -k , 6 + 2it[, respectively. Then Te{t) := n^(t) — Iig{t). 
Note that the matrix for Tg{t) in the basis of eigenvectors for Li{t) is diagonal with entries ±1, so detTg{t) = ±1 in- 
dependent of t. In the present case Li{t) = ^74 and it is easy to see that Tg = ±74 with sign "-|-" for — 7r/2 < 9 < Tr/2 
and "— " for w/2 < 6 < 37r/2. Now, inserting the above expressions for Li and Lo{t) into U{Li,Lo) , Sg{L-i,Lo) , and 
R{Li,Lo) , and inserting these into (6), and noting the following easy consequences of the definition (5) of Ha{t) 
(which follow from the vanishing of the traces of 71^ , 747,^ and the anti-hermitian matrix A^ipc, r)): 

Tr{HA{T) - fi) = -fiN (7) 
TriTeiHAiT) - ^)) = Tr{±^^{HA{T) - /x)) = ±mN (8) 

we arrive after a little calculation at the result 

detDAin) = C e-^'*^/2 ^^^^^ ^ i(A{f3)) (9) 

where 

UAiP)^Te'-fo"^^^^'^'' (10) 

and C = e^l^"^^/"^. (T=T-ordering and we have used the fact that A'' = dimW is formally even, since spinor fields have 
an even number of components, hence (—1)^ = 1-) We will discuss below how variants of this approach can produce 
different expressions for the overall factor C in (9). In the present case C = 6=*=^'"^/^ is gauge field-independent - 
its only effect is to give an overall shift in quantities obtained by taking derivatives of logZ with respect to /3 or m 
- so the the factors C e~^'^^/^ in (9) are physically inconsequential and can be absorbed into a normalisation of the 
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QCD partition function; hence the physics of the fermion determinant is described solely by det{l + e^'^L(A{(^))- Thus 
(9) is a dimensionally reduced expression for the fermion determinant, showing that the fermion action in (1) can be 
replaced by a non-local action for fermion fields living in one spacetime dimension lower (i.e. living only in the spacial 
volume). Moreover, it provides the promised generalisation of (3) to the case of arbitrary background gauge field: In 
the free field case, Uo{P) = e^^^° and it is easy to check that (9) reproduces (3) with Co = C. 

The result (9) also provides another way to relate the functonal integral and Hamiltonian frameworks for quantum 
field theories with fcrmions, alternative to the usual relation based on transfer matrices which requires a time discreti- 
sation [11]. In the free field case this is easily seen using a standard algebraic rewrite of the determinant in the right- 
hand side of (9): Defining the extension of an operator Q on to an operator on the Clifford algebra ®^o^^^ 
(~fermionic Fock space) by Q^(V-'iA- • -AV-p) = (QV'i)A- • •A(Qi/'p) we have det(l+e-''(^«-'')) = TrFock{{e~'^''"°~^^T) 
and the latter can be identified with Tr{e~^^^~'^'^^) where 7i and J\f are the fermionic Hamiltonian and number op- 
erators, respectively, on the Fock space. (In this description the vacumm state has unfilled Dirac sea, i.e. it consists 
of no particles of either positive or negative energy. Vacimm subtractions then need to be done "by hand" to obtain 
physical results.) Relating the functional integral and Hamiltonian approaches via (9) in the general case where gauge 
fields are present is more involved and will be discussed in a seperate paper. 

The above derivation of (9) is formal since W = {V-'(x)} is infinite-dimensional. However, a rigorous version can be 
obtained by putting the spacial volume on a lattice, i.e. by working in a continuous time — lattice space setting where 
W is the finite-dimensional vectorspace of spinor fields living on the sites of the spacial lattice. The Dirac operator in 
that setting was discussed in [5] (see eq.(7) of that paper); it is obtained from the above expression (4) (5) for -Da(/j,) 
by replacing dk + (=the covariant finite difference operator constructed with the link variables of the gauge 

field on the spacial lattice) and m — > Ma{t) = m + ^A^^gg(r) , where j^^fpacei"^) is the spacial Wilson term as 
defined in [5]. (Wc contirnic to denote the gauge field by A in the continuous time lattice space setting even though 
its spacial components are link variables.) The derivation of the dimensionally reduced expression for the fermion 
determinant then goes through as above but with one small change: m — > M^(t) in the right-hand side of (8), and 
consequently C = 6=*=^™^/^ gets replaced in (9) by 

C±(^)=e^*r'^'-^-M<^- (11) 

This modification is significant: C±{A) is gauge field-dependent and the indeterminacy of the sign in the argument 
of the exponential therefore constitutes an inconsistency in continuous time — lattice space QCD when the fermion 
determinant is defined via zeta-regularisation. There is a way to avoid this inconsistency though. Since "0 — "0*74 ths 
fermion determinant can be expressed as det{j4DA{^i)) (in fact this is the starting point for the usual evaluation of 
the fermion determinant via Matsubara frequency summation in the free field case). Although formally det{'^iDA{lj)) 
coincides with detDA{lj), the expressions obtained from zeta-regularisation turn out to be different. The difference 
originates from a difference in Tg : For 74Z)a(m) we find Tg = ±1 (with the sign depending on 9 in the same way as 
previously), and consequently eq.(8) becomes Tr{T0{HA{T) — fi)) = Tr{±{HA{T) — = '^\iN . The result is that 
det{^I^D Ai,^^)) is given by (9) with C = e^^^^^"^ , both in the formal continuum setting and in the regularised (spacial 
lattice) setting. Since this C is gauge field-independent the sign indeterminacy does not matter and the factor can be 
absorbed into a normalisation of the QCD partition function. 

Thus in the continuous time — lattice space setting the overall factor C in the dimensionally reduced expression (9) 
for the fermion determinant depends on choices made in the evaluation: it can be either C+{A) , C-{A) or e^^i^^l"^ . 
Clearly, the question of whether these choices are physically equivalent corresponds to the question of whether or 
not the gauge field-dependent factor C±(j4) is physically consequential. This same question has recently arisen in an 
investigation of universality in Lattice QCD [5]. In fact (9) has a very similar structure to a Lattice QCD fermion 
determinant expression obtained previously by an algebraic method in [12] (see also [13]). Its continuous time limit 
was evaluated in [5] for various versions of the Wilson/naive lattice fermion formulations, and, depending on the choice 
of lattice formulation, the limit was found to coincide with (9) with C = C-|- (^) or C = 1 (modulo some physically 
inconsequential factors). The factor C+(yl) was thus seen to represent a "universality anomaly" in that case. 

How serious a problem is this in the present case? Zeta-regularisation has become a rather standard technique 
which seems to have always given sensible, consistent results in the cases where it has been applied (see, e.g., [14]), 
while lattice regularisation is a well-established and successful approach at least in QCD. It would therefore be an 
unpleasant shock to find that the combination of these regularisations in the present setting leads to a physically 
ambiguous or inconsistent result for the fermion determinant. To determine whether this is actually the case one 
needs to determine the physical significance, or lack thereof, of the troublesome factor C±{^A). To this end, recall that 
^a{j) = m -|- ^^%acjy) ill (11)- One could argue that this factor is physically inconsequential when one goes on 
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to take the continuous space limit (a' ^ 0) since the spacial Wilson term ^A^^j^,. formally vanishes in this limit. 
This is a delicate issue though, since Tr^Agpace actually diverges in this limit (the largest eigenvalue is ~ ^). It is 
tempting to interpret this divergence as being due to the spacial "fermion doubler" modes, which get masses ~ ^ 
from the spacial Wilson term and are supposed to decouple in the spacial continuum limit. But the situation may be 
more complicated than this and further study is required to clarify this issue. 

In any case, the dimensionally reduced expression (9) sheds light on the /3- and /z-dependence of the fermion 
determinant and wc study this in the remainder of the paper, focusing in particular on the /i-dependence. Further 
details of our calculations, along with some additional results, are given in [15]. We begin by noting some general 
features. The operator Ua{0) given in (10) has the properties 

m*)-^ = (7475)-' U{I3) (7475) (12) 

detUAiP) = 1 (13) 

The derivation of the first of these is postponed to [15]; it immediately gives |dei Z//a(/3)P = 1, but to de- 
rive (13) we need to show the absence of a phase factor and this requires a direct calculation: -j^ detU{l3) = 

detU{l3)Tr{U{l3)-^ -^U{13)) = detU{^)TT{U{^)-\-H{^))U{P)) = -detU{P)Tr H{^) ; this vanishes by (7), hence 
detU{f3) = detU{0) = 1 as claimed. Using (12)-(13) it is straightforward to see that the dimensionally reduced 
expression (9) satisfies the standard relation 

detDAinT = detDA{-n*) , (14) 

(which follows at the formal level from Z)a(/u)* = 75 Da{—^i*) 75), showing that detDA{^l) is complex-valued in general 
but real for purely imaginary /z. 

Specialising to real we consider now the large /x limit of the fermion determinant. In this limit the fermion 
lagrangian in (1) is dominated by the term fi'ip"/4^'>p ; thus, formally, the dynamical fermion effects disappear from 
the QCD partition function in this limit and the thermodynamics is that of the pure gauge theory. This should 
correspond to the fermion determinant becoming independent of the gauge field in the large /j, limit, and the question 
of whether this happens can be investigated directly from the dimensionally reduced expression derived in this paper: 
As a consequence of (13) we have det{l + eJ^^^ Ua{i3}) « e^^^^ in the large fi limit, and it follows that in the formal 
continuum setting the expression (9) does indeed become independent of the gauge field. However, in the regularised 
(spacial lattice) setting a residual gauge field dependence remains, contained in the notorious factor C±{A). Thus the 
intuitive expectation of gauge field independence in the large fi limit is not entirely realised in this case. 

To study the /z-dependence of the fermion determinant in more detail it is useful to re-express (9) in terms of the 
eigenvalues of 74.0^(0) = + Ha{t) as follows. Consider the eigenvalue equation 

^ + HA{r)nr) = {X+jmr) , *(/?) = -*(0) (15) 

A little analysis shows that the solutions are given by *(r) = e*^^+^^'^ ZiA(r) *(0) with Z//a(/3) *(0) = e~'^^*(0). 
Prom this we see that the eigenvalues of ^ -|- Ha{t) come in equivalence classes {{A -|- ^(2n -|- l)}nez} which are 
in one-to-one correspondence with the eigenvalues {e"'^'^} of Ua{0)- The representative A for the equivalence class 
{A -t- ^(2n + l)}nez can be fixed by the condition Im{X) G [— 7r//3, 7r//3[. A known property of the eigenvalues that 
we exploit in the following is that the A's with Re{\) ^ come in pairs (A, —A*) (a derivation of this is provided in 
[15]). Using this in (9) we obtain 

detDAifJ.) =CY[ e-^^'^ (1 + e'^" e"'^^) 

A 

= C[ n e-P^/^l+eP^^-^^)\[ W 6-/^^(1 +e'^('^-^))(l + e'^(''+^*))] (16) 

Re{X)=0 Re(\)>Q 

Note that the last square-bracketed factor can be rewritten as ni?e(A)>o + e~''('^~''^)(l -|- e~^'^^'^'^^), which 

has the same structure as the earlier free field expression (3). This can now be used to obtain an expression for the 
fermion determinant in the limit of low temperature (i.e. large /?). After some manipulations we find 

det£>A(M)''«"Ce-'5^-(^)e(^) [] e'^d'^l-I^^^D/^ e-'^'w (^^0) (17) 

\Re(\)\<M 

detDA{0)'^^°° Ce-^^^"^^h{A) J| 2cos{P Im{X)/2) (18) 

Re{\)=0 
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where EgeaiA) X]fie(A)<o -^'^(A) is the sum of the negative quasi-encrgies (this is the natural generalisation of the 
sum of energies of the negative energy states in the Dirac sea arising in the free fermion case), and 

e{A):=l[e-'P'"'^^'>/\ (19) 

A 

The latter is a gauge field-dependent sign factor: e{Af = f]^ e-*^^'"^^) = Ha ^"^^ = detUA{P) = 1 
by (13). We now consider this sign factor in more detail. Setting = in (16) we find detDA{0) = 
e(^) nfle(A)=o2cos(/3/m(A)/2) nfle(A)>o + e-'^^)\'^ el^ which shows that e{A) determines the sign of 
detDAiO). (Note that the condition Im{X) G [— tt//?, n/ 13[ implies cos{(3 Im{X)/2) > 0. If wc replace A ^ A + i27rn//3 
then for i?e(A) = the change in sign of e(A) is compensated by a change in sign of cos(/3/m(A)/2), while for Re{X) > 
the replacement leaves e{A) unchanged since A enters in e{A) through the factors e"'''^"*^'^)/^ e"*''^™^"'^ = 
g-i;3 7m(A) -J Formally, the sign of detDA{0) is given by the sign of the product of the real eigenvalues of Da{0); 
this follows from the formal relation det{DA{0) - z*) = det{DA{0)* - z)* = det{j5DA{0h5 - z)* = det{DA{0) - z)* 
which shows that the non-real eigenvalues come in complex conjugate pairs. In particular, if the real eigenvalues 
arc > then detDA{0) > 0. In the present regularised (lattice space) setting the spacial Wilson term is a positive 
operator; consequently, assuming m > 0, the real parts of the eigenvalues of Da{0) are positive; hence one would 
expect detDA{0) > and therefore e(^) = 1 for all gauge fields A. We now show rigorously that this is indeed the 
case. It is a nontrivial problem to show this directly from the definition (19); however, it can be inferred indirectly 
as follows. Modulo some (positive) inconsequential factors, the dimensionally reduced expression (9) was found in 
Eq.(20) of [5] to arise as the continuous time limit limo_>o a^^^ detDa^ where Da^ (with a = ii + in/P) is a lattice 
time — lattice space version of D{fi). Now, since D^^ is a finite-dimensional matrix and satisfies (Da^)* = 75 Da^ 75 
for fi = (i.e. a = ztt//?), the preceding reasoning rigorously implies detDa^ > 0. Its a ^ limit must therefore 
also be positive, hence detDA{0) > and therefore e(^) = 1 as claimed. (We remark though that the sign can be 
negative when m < 0, since in this case negative real eigenvalues are possible. This situation occurs in Lattice QCD 
with Wilson fermions: to study the chiral limit m is tuned to negative values and there are so-called "exceptional 
configurations" for which the Dirac operator then has negative real eigenvalues.) 

The large (3 expressions (17)-(18) have a remarkably simple form. Besides the natural factor e~^^""^, and the sign 
factor e(A), which is positive for all A when m > but which can be negative in certain gauge backgrounds when 
m < 0, the /z = expression (18) consists only of a simple factor associated with the purely imaginary A's, while in 
the /U 7^ case (17) there is a further simple factor associated with the A's for which |iie(A)| < In particular, we 
see the following: The fermion determinant at large /3 is independent of /j. for < ndA) where 

fio{A) :=min{\Re{X)\}. (20) 

This property was also found in [4] from an expression for the ratio detD a{ij) / detD a{Q) derived via the Matsubara 
summation method. By considering (a suitable integral over) the pion propagator, miri{i.io{A)} can be identified 
with TO7r/2 [4] (where the minimum is taken over the statistically significant gauge fields; similar considerations are 
standard in the Lattice QCD setting; sec, e.g., [12].) It follows that the zero temperature QCD partition function 
is independent of /i for jjl < m7r/2, as discussed at the beginning of this paper. The results (17) (18) reproduce the 
determinant ratio expression Eq. (6) of [4] (where it was implicitly assumed that i?e(A) ^ for all A), and extend 
that result to the fermion determinant itself, revealing the presence of gauge field-dependent factors Q-fE^^aiA) ^^-^^ 
which cancel out in the determinant ratio. Having obtained the large (3 expression for the fermion determinant itself 
we are now in a position to investigate the subtle cancellation in the functional integral over gauge fields in (1), 
which, as discussed earlier, is required to explain the ^^-independence the zero-temperature QCD partition function 
for m7r/2 < /i < m^/S. We hope to gain more insight into this in future work via the program outlined earlier. 

Besides the aforementioned program for gaining insight into the QCD phase transition at zero temperature and 
finite density there are other places where the results obtained here may be useful. For example, it would be interesting 
to see, if possible, how the high density effective theory (HDET) of Hong and Hsu [16], where the effective fermion 
determinant is positive, emerges from the dimensionally reduced determinant expression (9). This may shed light 
on possible limitations of that approach. We conclude by mentioning another setting where the method introduced 
in this paper for obtaining dimensionally reduced determinant expressions via a (partial) zeta-regularisation, using 
the formula Theorem 1 of [8], may also have interesting applications. Following the Randall-Sundrum proposal [17] 
there has been much interest in gauge theory models where the spacetime includes a fifth dimension interval, with the 
geometry of the interval being warped so that the spacetime is a slice of AdS^. In particular, a model has recently 
been proposed in this framework where the fermion masses arise not from the condensation of a Higgs field (the 
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standard mechanism) but as a consequence of the boundary conditions on the fifth dimension interval [18]. It would 
be very interesting to derive a dimensionally reduced expression for the fermion determinant in this model - this 
would reveal a non-local 4-dimensional fermion theory to which the 5-dimensional model is equivalent. (This proposal 
is reminiscent of the derivation of the overlap Dirac operator in lattice gauge theory [19] from a dimensional reduction 
of the fermion determinant of a 5-dimensional domain wall model [20] . Since the overlap Dirac operator has proved to 
be very interesting and useful both conceptually and for practical purposes, one may expect the same for the non-local 
Dirac operator in 4 dimensions obtained from a dimensionally reduced expression for the fermion determinant of the 
5-dimensional model of [18].) The boundary conditions there may be too complicated for a simple application of the 
determinant formula of [8] such as the one in the present paper though; the result of [8] may first need to be extended 
to more general boundary conditions. (We remark that more general boundary conditions have been considered in 
[9], although the results there are mostly for determinant ratios.) The effect of the warped geometry of the AdS^ slice 
on the fermion determinant will also need to be dealt with, although this appears not so difficult. We also remark 
that in light of the AdS/CFT correspondence [21] the dimensionally reduced expression for the fermion determinant 
in this model can be expected to contain CFT structure and it would be interesting to uncover this. 
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